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Abstract.  Dynamical  theories  of  the  Earth’s  rotation  like  SMART97 
(Bretagnon  et  al.,  1998)  are  to  be  considered  in  a  DGRS  (dynamically 
nonrotating  geocentric  reference  system)  (Brumberg  et  al.,  1996).  Such  a 
theory  gives  the  explicit  expressions  in  terms  of  TCG  (Geocentric  Coordi¬ 
nate  Time)  of  three  Euler  angles  relating  a  DGRS  to  the  ITRS  (Interna¬ 
tional  Terrestrial  Reference  System).  These  angular  quantities  together 
with  their  TCG  derivatives  enable  one  to  get  all  Earth’s  rotation  param¬ 
eters.  At  the  same  time,  the  analysis  of  observations  result  in  the  val¬ 
ues  for  slightly  different  angles  and  their  TCG  derivatives  characterizing 
the  relationship  between  the  ITRS  and  a  KGRS  (kinematically  nonro¬ 
tating  geocentric  reference  system).  The  differences  between  these  two 
sets  of  six  quantities  represent  kinematical  relativistic  corrections  (due 
to  geodesic  precession,  geodesic  nutation  and  luni-planetary  terms).  The 
paper  presents  these  differences  computed  by  means  of  the  VSOP87  series 
(Bretagnon  and  Francou,  1988).  In  particular,  in  analysing  observations 
at  the  microarcsecond  level  these  expressions  will  permit  an  experimen¬ 
tal  check  of  geodesic  precession  in  a  more  direct  manner  than  it  is  done 
nowadays  (Bertotti  et  al.,  1987). 


1.  Introduction 

The  aim  of  this  paper  is  to  supply  numerical  values  for  the  quantities  involved 
in  the  relativistic  barycenric  and  geocentric  reference  systems  (RSs)  defined  in 
Brumberg  et  al.(  1996)  and  Brumberg  (1997a).  Using  the  VSOP87  series  for 
the  motion  of  the  major  planets  (Bretagnon  and  Francou,  1988)  this  work  was 
started  in  Brumberg  et  al.(  1992)  by  computing  the  expressions  for  the  geodesic 
rotation  vector  and  its  derivative  and  in  Bretagnon  et  al.(  1997)  by  computing 
the  differences  of  the  Euler  angles  for  the  geocentric  rotation  matrix  relating 
the  dynamical  and  kinematical  systems.  Our  aim  is  to  compute  the  differences 
of  the  Euler  angles  derivatives  as  well.  In  anticipating  the  space  astrometry 
projects  of  microarcsecond  accuracy  it  enables  one  to  consider  completely  the 
effect  of  the  geodesic  rotation  on  the  Earth’s  rotation  parameters. 
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2.  Barycentric  and  geocentric  RSs 

In  dealing  with  RSs  we  use  the  same  notation  as  in  Brumberg  et  al.(  1996)  and 
Brumberg  (1997a),  i.e.  B  -  barycentric,  G  -  geocentric,  V  -  VLBI,  C  -  ecliptical, 
Q  -  equatorial,  D  -  dynamical,  K  -  kinematical,  +  -  rotating.  At  the  barycentric 
level  we  have  three  systems  BRSV,  BRSC  and  BRSQ  with  timescale  t=TCB  and 
spatial  coordinates  x  =  (ad),  xc  =  (%c)  and  xq  =  (xq),  respectively.  At  the 
geocentric  level  each  of  these  systems  generates  by  means  of  the  BRS^GRS  two 
geocentric  RSs,  dynamically  or  kinematically  nonrotating  ones  with  respect  to 
the  generating  BRS.  As  a  result,  there  will  be  six  geocentric  RSs,  i.e.  GRSV, 
GRSC  and  GRSQ  for  D  and  K  versions  with  timescale  n=TCG  and  spatial 
coordinates  w  =  (io!),  wq  =  (wq)  and  wq  =  ( Wq ),  respectively.  If  necessary,  we 
will  distinguish  D  and  K  versions  by  writing  i»'  with  q  =  1  for  version  D  and  q  = 

<7 

0  for  version  K.  One  more  geocentric  system  is  GRS+  rotating  with  the  Earth  and 
having  spatial  coordinates  y  =  (yr).  By  identifying  BRSV  and  GRS+  with  ICRS 
and  ITRS,  respectively,  the  problem  is  to  determine  all  relationships  between 
these  systems  enabling  one  to  obtain  an  unambiguous  relativistic  interpretation 
of  Earth’s  rotation  parameters  as  well  as  related  astronomical  concepts. 

At  the  barycentric  level  all  three  RSs  are  mutually  related  by  the  constant 
rotation  matrices 

xc  =  Pcx  ,  XQ  =  PQx  .  (2.1) 

Very  approximately 

(10  0  \ 

Pq  =  E,  Pc  =  I  0  cose  sme  I  ,  (2.2) 

\  0  —  sin  e  cos  e  / 

where  E  stands  for  the  unit  matrix  and  e  is  the  mean  obliquity.  This  constant 
rotation  is  conserved  at  the  geocentric  level  in  the  relationships  between  V,  C 
and  Q  GRSs  of  the  same  type  (with  respect  to  D  or  K  versions).  The  matrix 
Pc  may  be  determined  from  the  comparison  of  the  VSOP87  planetary  theories 
(Bretagnon  and  Francou,  1988)  with  observations  of  the  major  planets  within 
BRSV  (ICRS)  background.  Matrix  Pq  is  to  be  determined  at  the  geocentric 
level  by  comparison  of  the  Earth’s  rotation  theory  like  SMART97  (Bretagnon 
et  al. ,  1998)  with  observations  within  GRS+  (ITRS)  background  (see  below). 
The  relationship  with  GRS+  involves  the  Earth’s  rotation  matrix  P(u) 

y  =  p(u)  WC  =  P(u)PC  W  =  P(u)PCPq  Wq  .  (2.3) 

q  q  q  q  q  ^  q 

The  Earth’s  rotation  matrix  may  be  expressed  in  terms  of  three  Euler  angles  to 
be  determined  by  some  Earth  rotation  theory.  Following  the  traditional  choice 
by  Tisserand  (1891)  and  Smart  (1953)  we  will  use  here 

P(u)  =  D3(¥))jD1(-0)jD3(-^),  P(u)  =  D^^Dri-e^D^-^K),  (2.4) 

1  0 

with  elementary  rotations 


/I 

0 

.°  \ 

(  cos  (3 

0 

—  sin  (3 

Di(a)  =  0 

cos  a 

sin  a  , 

D2(P)  = 

,° 

1 

0 

\o 

—  sin  a 
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\  sin  [3 

0 

cos  (3 

Relativistic  Corrections  for  ERPs 


295 


(cos  7  sin  7  0  \ 

—  sin 7  cos  7  0  )  .  (2-5) 

0  0  l) 

The  relationships  between  the  D  and  K  systems  involve  the  geodesic  rotation 
matrix  F 

w  =  (E-c~2F)w  ,  wc  =  (E  -  c~2Fc)wg  ,  wQ  =  (E  -  c~2Fq)wq,  (2.6) 
and 

P{u)  =  P{u)(E  -  c~2Fc),  (2.7) 

1  0 

with 

FC  =  PcFPl ,  Fq  =  PqFPq  .  (2.8) 

Instead  of  rotation  matrix  F  one  may  consider  its  equivalent  vector  representa¬ 
tion 

F13  =  £lJkFk  ,  £ijk  =  \{i  -  j){j  -  k){k  -  i),  (2.9) 

so  that  the  matrix  product  in  (2.6)  may  be  reduced  to 

Fw  =  —SijkF3  wk  .  (2.10) 

The  components  of  Fk<  and  Fk<  have  been  evaluated  in  Brumberg  et  al.(  1992). 
The  relationships  between  barycentric  and  geocentric  RSs  are  based  on  the 


BRS- 

-^GRS  transformation  (Brumberg,  1995  and  references  therein) 

u  =  t  —  c  2[A(t)  A  vete\  T  •  •  •  5 

(2.11) 

wi 

=  rlE+  c"2{[|44  +  +  Dlk(t)]rkE  +  Dlkm (t)rkEr%}  +  . . . , 

■  (2-12) 

with 

rlE  =  xl  -  xlE(t),  vlE  =  ilE(t), 

(2.13) 

A(t)  =  \v2e  +  Ue(*e), 

(2.14) 

and 

Dik(t)  =  StkUE(xE).  Dl3k(t)  =  \(6l3akE  +  6tka3E  -  FjkFE ), 

(2.15) 

where  Ue(x-)  stands  for  the  Newtonian  potential  of  all  solar  system  bodies  except 
the  Earth,  xlE,  vE  and  aE  being  Earth’s  BRS  position,  velocity  and  acceleration, 
respectively.  This  transformation  is  written  for  V  versions  of  RSs.  For  C  or  Q 
versions  one  should  convert  all  V  quantities  to  C  or  Q  quantities,  respectively. 

Let  t*  be  the  BRS  moment  of  time  corresponding  to  an  event  with  the  GRS 
coordinates  (u,wl  =  0)  (Klioner  and  Voinov,  1993).  Then  u  and  t*  are  related 
by  the  time  equation 

u  =  t*  -  c~2A(t*)  + ...  .  (2.16) 


=  c~2vkEwk 


Hence, 


t  —  t 


(2.17) 
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Expanding  the  right-hand  member  of  (2.12)  in  the  vicinity  of  t*  one  gets  the 
inverse  GRS^BRS  transformation 

t  =  u  +  c  2[H(m)  +  zEwk]  +  . . . ,  (2.18) 

X i  =  wi  +  4M  +  c“2  [(|44  -  (lplk  -  Dlk)wk  -  Dtkmwkwm ]  +  .  . .  ,  (2.19) 

with  functions  zE  =  zE(u )  characterizing  the  motion  of  the  geocenter  in  terms 
of  TCG  and  determined  by 


4M  =  4(0-  (2-2o) 

This  inverse  transformation  is  used,  for  example,  to  transform  GRS  space-time 
coordinates  of  terrestrial  ground  stations  into  their  BRS  space-time  coordinates. 


3.  K— D  for  Euler  angles 


The  first  of  relations  (2.3)  for  the  dynamical  version  q  =  1  represents  a  relation¬ 
ship  between  GRS+  (ITRS)  and  DGRSC.  The  Euler  angles  g>,  ip  and  6  together 
with  their  TCG  derivatives  result  from  the  solution  of  the  Earth’s  rotation  equa¬ 
tions  written  in  the  DGRSC.  This  solution  may  be  thought  of  as  consisting  of 
two  parts,  a  Newtonian  one  corresponding  to  the  formal  Newtonian  equations 
of  the  Earth’s  rotation,  and  relativistic  correction  terms  generated  by  the  rela¬ 
tivistic  contributions  to  the  right-hand  sides  of  the  Earth’s  rotation  equations. 
These  dynamical  corrections  8(p,  Sip  and  SO  as  well  as  their  TCG  derivatives  are 
neglected  here  due  to  their  small  size  in  the  DGRSC  (Brumberg,  1997b,  Klioner, 
1997,  and  references  therein).  For  the  Newtonian  solution  of  the  Earth’s  rotation 
problem  we  choose  here  the  SMART97  theory  (Bretagnon  et  al. ,  1998).  On  the 
other  hand  astronomical  observations  of  the  Earth’s  rotation  may  be  interpreted 
as  performed  in  KGRSC  resulting  in  the  angles  g>K,  ipK  and  Ok  and  their  TCG 
derivatives.  They  differ  from  their  dynamical  counterparts  due  to  the  geodesic 
rotation  vector  FK .  Using  (2.4)  and  (2.7)  it  follows  that 

_2 

<P~<PK  =  --z—z  (Fq  s’mip  +  Fq  cosip)  ,  (3.1) 

sin  0  v  ' 


and 


0  —  Ok  =  c  2  ^ Fq  cos  ip  —  Fq  sin  ip^j  , 


4>  ~  tpK  =  C 


-2 


o  cos  0  ,  o 

4)  -  {f};  sin  Ip  +  Fl  cos  ip 


sin  0 


(3.2) 

(3.3) 


In  this  way  the  matrix  P(u)  is  known  (for  both  D  and  K  versions)  from  the 
Earth  rotation  theory.  Having  found  the  constant  matrix  Pq  from  the  planetary 
motions  one  may  determined  the  relationship  between  GRS*  and  GRSV  from 
the  second  of  relations  (2.3).  The  main  plane  of  GRS+  is  the  equator  of  date. 
Having  determined  the  mutual  orientation  of  GRS+  and  GRSV  one  may  find 
the  constant  matrix  Pq  defining  the  main  plane  of  GRSQ  as  a  fixed  plane  for 
some  definite  epoch.  The  inverse  transformation  GRSQ^BRSQ  enables  one  to 
fix  the  BRSQ. 
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Expressions  (3.1)-(3.3)  were  evaluated  first  in  (Bretagnon  et  al. ,  1997)  using 
the  series  of  SMART97.  The  three  main  RSs  considered  in  (Bretagnon  et  al., 
1997)  are 

Oxyz  =  GRSC,  O^C  =  GRS+,  Ofi)C  =  D3(a)GRS+ , 

a  being  the  longitude  of  the  axis  associated  with  the  principal  moment  of  inertia 
of  the  Earth.  This  means  that  the  system  0 is  a  terrestrial  system  with 
principal  axes  of  inertia.  Therefore,  Bretagnon  et  al.(  1997)  angles  u  and  ip 
are,  respectively,  —  9  and  —  ip  of  Brumberg  et  al.(  1996)  and  Brumberg  (1997a) 
whereas  angle  ip  has  the  same  meaning  in  these  papers.  The  expression  for 
ipK  ~  A  given  in  Bretagnon  et  al.(  1997)  should  be  taken  with  the  opposite  sign. 


4.  K— D  for  Euler  angles  derivatives 

To  evaluate  the  difference  in  the  values  of  the  Earth’s  rotation  parameters  (ERP) 
in  dynamical  and  kinematical  RS  one  should  get  along  with  ip  —  ipK,  ip  —  ipK  and 
9  —  9 k  their  derivatives  <p>  —  CpK,  'ip  —  *Pk  and  9  —  9k-  They  may  be  computed 
by  means  of  formulas 

_2 

¥>  -  ipK  =  \tc  sin  +  Fc  cos  +  (0  -  Ok)'>P  +  (<P~  <PI<)9  cos  9] ,  (4.1) 

sin  9  L  J 

•  •  •  c~2  r  •  -l 

ip  —  ipK  =  c  2Fq  — : — -  \  (Fq  sin  ip  +  F2  cos  ip)  cos  9  +  (9  —  9k  )'>P  cos  9  +  (<p>  —  <pk)0 J  ? 

(4.2) 

and 

9  —  9k  =  c~2[Fq  cos  ip  —  F2  sin  ip  +  (ip  —  <pK)i>  sin  8],  (4-3) 

using  (3.1)-(3.3)  and  SMART97  theory  (Bretagnon  et  al.,  1998)  for  the  angles 
ip,  9  and  their  derivatives.  The  initial  terms  of  the  series  for  (4.1)-(4.3)  are  given 
in  the  Appendix. 


5.  Earth’s  rotation  parameters 

Representation  (2.4)  means  splitting  P(u)  into  matrix  S  of  the  diurnal  rotation 
and  matrix  N  of  the  precession  and  nutation,  i.e. 

P(u)  =  SN ,  S  =  D3(<p),  N  =  D1(-9)D3(-iP)  (4.1) 

(ignoring  for  simplicity  the  motion  of  the  Earth’s  poles).  This  splitting  involves 
the  representation  of  the  total  Earth’s  rotation  velocity  w'  (GRSC  components 
of  the  angular  velocity  of  rotation  of  GRS+  with  respect  to  GRSC)  in  two  parts 

uf1  =  uF-jk  +  Njiojg,  (4-2) 


with 


(° 

^5=  0 

\<P 


(—9  cos  ip 
9  sin  ip 
-iP 


(4.3) 
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By  applying  these  relations  to  the  D  and  K  versions  one  gets 


us  -  u  s 


0 

0 

A  -  (pK 


(4.4) 


and 


UN  ~  UN 
1  0 


^  (ip  —  iPk)0  sin  ip  —  {0  —  Or)  cos  ip 

(ip  —  iPk)0  cos  ip  +  (0  —  Ok)  sin  ip 

V  ~{i>  ~  P’k) 


(4.5) 


Relation  (4.4)  characterizes  the  difference  of  the  determination  of  the  sidereal 
time  in  dynamical  and  kinematical  systems.  Relation  (4.5)  gives  the  difference  of 
precession-nutation  contributions  in  the  Earth’s  angular  velocity  in  dynamical 
and  kinematical  systems.  The  initial  terms  of  the  series  representing  the  two 
first  rows  of  matrix  (4.5)  are  reproduced  in  the  Appendix. 


6.  Conclusion 

Using  VSOP87  planetary  theories  and  SMART97  Earth  rotation  theory  we  have 
computed  the  relativistic  kinematical  differences  in  the  Euler  angles  and  their 
TCG  derivatives  relating  the  ITRS  and  the  geocentric  ecliptical  reference  system 
GRSC  in  dynamical  and  kinematical  versions.  For  future  observations  at  the 
microarcsecond  level  it  enables  one  to  consider  more  rigorously  the  Earth’s  ro¬ 
tation  parameters  in  different  reference  systems.  For  the  present  analysis  it  was 
sufficient  for  us  at  the  barycentric  level  to  relate  the  ICRS  with  equatorial  and 
ecliptical  reference  systems  just  by  means  of  simple  matrices  (2.2).  The  initial 
terms  of  the  series  determined  here  are  given  in  the  Appendix.  The  complete 
series  are  available  by  request  to  the  second  author. 

In  view  of  the  present  analysis  it  is  interesting  to  consider  the  possibility 
of  checking  geodesic  precession  directly  by  comparing  ERP  referred  to  DGRSC 
and  KGRSC.  Geodesic  precession  has  been  tested  only  implicitly  by  analyzing 
the  lunar  perigee  motion  with  the  aid  of  LLR  and  VLBI  observations  (Bertotti 
et  cil. ,  1987). 

Acknowledgments.  The  first  author  appreciates  the  support  of  the  Rus¬ 
sian  Foundation  of  Fundamental  Researches  (Grant  No.  99-02-16793). 
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Appendix 

We  reproduce  below  the  initial  terms  of  the  series  for  p>  —  p>K,  —  fin,  0  —  9k, 
fi  —  ipK ,  —  6  —  &K  and  first  two  rows  of  matrix  (4.5).  The  time  t  is  measured 

in  Julian  years  from  J2000.  These  series  involve  12  angular  variables,  i.e.  8  mean 
planetary  longitudes  with  respect  to  the  fixed  equinox  J2000  A;  (1  <  i  <  8),  the 
Delaunay  arguments  D,  F,  l  of  the  Moon  and  the  mean  angle  of  rotation  of 
the  Earth  <p.  These  arguments  expressed  in  radians  are  represented  by  linear 
functions  of  time  as  follows: 

Ai  =  4.40260867435  +  26087.90314157421 , 

A2  =  3.17614652884+  10213.28554621101 , 

A3  =  1.75347029148  +  6283.07585114551 , 

A4  =  6.20347594486  +  3340.61242669981 , 

A5  =  0.59954632934  +  529.69096509461 , 

A6  =  0.87401658845  +  213.29909543801 , 

A7  =  5.48129370354  +  74.78159856731 , 

A8  =  5.31188611871  +  38.13303563781 , 

D  =  5.19846640063  +  77713.77144818041 , 

F  =  1.62790513602  +  84334.66157178371 , 

l  =  2.35555563875  +  83286.91424771471 , 

Cp  =  4.89496121282  +  2301216.75365153651 . 

These  values  are  slightly  different  from  those  given  in  (Bretagnon  et  al. ,  1998). 
The  following  three  series  representing  the  differences  p>  —  p>K,  —  fin,  9  —  9k 
are  reproduced  here  with  some  more  terms  as  compared  with  (Bretagnon  et  al ., 
1997),  their  coefficients  being  expressed  in  pas: 

cp  -  (pK  =  -2.991-  54771.0312  -  802.0613  + 

+  3.28  sin(A3  +  D  -  F)  -  0.24  sin(2A5  -  5A6)  -  0.36  cos(2A5  -  5A6)  + 

+  0.07  sin(8A2  -  13A3)  -  0.04  cos(8A2  -  13A3)  +  0.02  sin(3A2  -  5A3)  + 
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+  0.07  cos(3A2  —  5A3)  —  0.05  sin(A3  +  F)  —  0.04  sin(A3  —  F)  + 

+  0.03  cos2A5  +  0.02  cos2A6  +  0.02  cos  A5  —  0.01  cos(A3  +  A5)  + 

+  7 [0 .20  sin  A3  +  0.85  cos  A3  +  0.80  cos(A3  +  D  —  F)  —  0.14  sin(2As  —  5Ag)  + 

+  0.05  cos(2A5  -  5A6)  -  0.02  sin(8A2  -  13A3)  -  0.02  cos(8A2  -  13A3)  + 

+  0.02  sin(3A2  -  5A3)  -  0.01  cos(A3  +  F)]  + 

+  l2  [—0.51  sin(A3  -\-  D  —  F  )  —  5.64  cos(A3  -\-  D  —  F  )  -\-  0.03  sin  2A3  +  0.34  cos  2A3  + 
+  0.05  sin  A3  —  0.02  cos  A3  +  0.06  cos(2A3  +  2D)  —  0.06  cos(2A3  +  2D  —  2 F)  + 

+  0.02  cos(2A5  —  5Ag)  +  0.01  sin  3A3  —  0.01  cos(A3  +  D  +  _F)]  + 

+  l3 [2.37 sin(A3  +  D  -  F)  -  0.20cos(A3  +  D  -  F)  -  0.24sin2A3  +  0.02cos2A3  - 

—  0.04  sin(2A3  +  2D)  +  0.04  sin(2A3  +  2D  —  2F)  —  0.02  cos  A3]  +  . . . ,  (A.l) 

ip  -  ipK  =  19198827.441  -  50386. 3212  -  754. 0913  - 

—  34.28  sin  A3  —  149.22  cos  A3  +  3.01  sin(A3  +  D  —  F)  —  1.73  sin  2A3  + 

+  0.84  cos  2A3  —  0.26  sin(2As  —  5Ag)  —  0.33  cos(2As  —  5Ag)  —  0.37  sin  D  — 

—  0.05  sin(4A3  —  8A4  -)-  3A3^  -)-  0.17  cos(4A3  —  8A4  -)-  3A5)  -)-  0.21  sin[A3  —  A5)  — 

—  0.13  sin(2A2  —  2A3)  +  0.09  sin(A2  —  A3)  +  0.09  sin(8A2  —  13A3)  + 

+  t[— 7.36  sin  A3  +  6.47  cos  A3  +  0.73  cos(A3  +  D  —  F)  +  0.19  sin  2A3  + 

+  0.15  cos2A3  —  0.13  sin(2A5  —  5Ag)  +  0.04  cos(2As  —  5Ag)  — 

—  0.02  sin(8A2  -  13A3)  -  0.03  cos(8A2  -  13A3)  +  0.02  sin(3A2  -  5A3)  - 

—  0.02  sin(4A3  —  8A4  +  3As)  —  0.01  cos(A3  +  _F)]  + 

+  l2  [—0.47  sin(A3  +  D  —  F)  —  6.15  cos(A3  +  D  —  F)  +  0.30  sin  A3  +  0.29  cos  A3  + 

+  0.04  sin2A3  +  0.35  cos  2A3  +  0.06  cos(2A3  +  2D)  —  0.06  cos(2A3  +  2D  —  2F)  + 

+  0.01  cos(2A5  —  5Ag)  +  0.01  sin3A3  —  0.01  cos(A3  +  D  +  +)]  + 

+  l3 [2.41  sin(A3  +  D  -  F)  -  0.20cos(A3  +  D  -  F)  -  0.25  sin  2A3  +  0.02  cos  2A3  - 

—  0. 04  sin(2A3  +  2D)  +  0. 04  sin(2A3  +  2D-2F)  -0.03  cos  A3]  +  . . . ,  (A. 2) 

9  -  0K  =  9.551  +  1954. Ill2  -  4721. 8013  - 

—  1.30  cos(A3  +  D-  F)  +  0.17  sin(2A5  -  5A6)  -  0.09  cos(2A5  -  5A6)  + 

+  0.02  sin(8A2  —  13A3)  +  0.03  cos(8A2  —  13A3)  —  0.03  sin(3A2  —  5A3)  + 

+  0.02  cos(A3  +  F)  +  0.02  cos(A3  —  F)  +  0.01  sin  2A5  + 

+  1[0.32  sin(A3  +  D  -  F)  -  0.03  sin(2A5  -  5A6)  -  0.04  cos(2A5  -  5A6)  - 

—  0.03  cos  A3]  + 

+  l2  [—1.82  sin(A3  +  D  —  F)  +  0.04  cos(A3  +  D  —  F)  +  0.14  sin  2A3  +  0.02  sin  A3  + 

+  0.08  cos  A3  +  0.02  sin(2A3  +  2D)  -  0.02  sin(2A3  +  2D-  2F )]  + 

+  13[— 0.04  sin(A3  +  D  —  F)  —  0.45  cos(A3  +  D  —  F)  +  0.07  cos  2A3  + 

+  0.01  cos(2A3  +  2D)  —  0.01  cos(2A3  +  2D  —  2F)\  +  . . . .  (A. 3) 

The  derivatives  ip  —  (px,  'ip  —  t/’/v,  9  —  9k  could  be  obtained  by  differentiating 
(A.1)-(A.3).  But  we  give  below  their  expressions  resulting  directly  from  (4.1)- 
(4.3).  Differentiation  of  (A.1)-(A.3)  was  used  only  for  checking  purposes.  The 
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series  for  the  differences  p  —  px,  P  ~  Pk-,  6  ~  @K  multiplied  by  106  are  as  follows: 

106(+  -  Pk)  =  —0.531  —  0.0112  +  0.0313  — 

—  0.02  cos(A3  +  F)  +  0.01  cos(A3  —  F)  — 

—  0.031  sin  A3 — 

—  t2  [0.05  sin(2A3  +  2D)  +  0.02  sin  2A3]  - 

—  t3  [0.03  cos(2 A3  +  2D)  +  0.01  cos  2 A3]  + 

+  0.0U4  sin(2A3  +  2D)  +  . . . ,  (A.4) 


106(>  -  pK)  =  93.08  -  0.491  -  0.0112  +  0.0213  + 

+  4.55  sin  A3  —  1.04  cos  A3  —  0.05  sin  2A3  —  0.11  cos2A3  —  0.14  cos  D  — 

—  0.02  cos(A3  +  F)  —  0.02  cos (D  +  /)  +  0.01  cos(A3  —  F)  + 

+  l[— 0.20  sin  A3  —  0.22  cos  A3  +  0.01  cos  2A3]  + 

+  12[— 0.05  sin(2A3  +  2D)  -  0.02  sin  2A3  +  0.01  sin(A3  +  D  +  F)  - 

-  0.01  sin(2A3  +  2 D  +  l)  -  0.01  sin(A3  +  D  -  F)] - 
—  I3 [0.03  cos(2A3  +  2D)  +  0.02  cos  2A3]  + 

+  0.0114  sin(2A3  +  2D)  +  . . . ,  (A. 5) 

1O6(0  -  0K)  =  0.021  -  0.0712  +  0.0212  cos(2A3  +  2D)+  ....  (A. 6) 

Finally,  we  reproduce  the  series  representing  the  first  two  rows  of  matrix  (4.5) 
multiplied  by  106: 

lO6^  -  =  -0.021  +  0.0712  + 

+12[— 1.76  sin(2A3  +  2 D)  -  0.02cos(2A3  +  2D)  -  0.77sin2A3  + 

+0.36  sin(A3  +  D  +  F)  —  0.35  sin(2A3  +  2D  +  /)  —  0.34  sin(A3  +  D  —  F)  + 

+0.07  sin(A3  +  +>  +  +  +  /)-  0.07  sin(2A3  +  4 D  -  l)  -  0.06  sin(2A3  +  4 D)  + 

+0.04  cos  3A3  -  0.05  sin(2A3  +  2D  +  21)  +  0.05  sin(2A3  +  2 D  -  l) - 

—0.03  sin(A3  +  D  —  F  +  /)  —  0.03  sin(A3  +  D  —  F  —  l)  —  0.01  sin(2A3  +  4 D  +  /)  + 

+0.01  sin(A3  +  3  D  +  +  —  /)  +  0.01  sin(2A3  +  /)  +  0.01  sin(A3  +  3  D  +  F1)]  + 

+13 [0.03  sin(2A3  +  2D)  —  0.86  cos(2A3  +  2D)  +  0.01  sin2A3  —  0.37  cos2A3  — 

-0.17  cos(2A3  +  2 D  +  l)  +  0.09  cos(A3  +  D  +  F)  -  0.08  cos(A3  +  D  -  F)  - 
-0.03  cos(2A3  +  4 D-l)-  0.03  cos(2A3  +  4 D)  -  0.02  cos(2A3  +  2D  +  21)  + 

+0.02  cos(2A3  +  2D  —  l)  —  0.02  sin3A3  +  0.02  cos(A3  +  D  +  F  +  l)]  + 

+14 [0.23  sin(2A3  +  2D)  +  0.02  cos(2A3  +  2D)  +  0.10  sin2A3  + 

+0.05  sin(2A3  +  2D  +  /)  —  0.01  sin(A3  +  D  +  F)  +  0.01  sin(A3  +  D  —  F1)]  + 

+15  [0.04  cos(2A3  +  2F>)  +  0.02  cos  2A3]  +  . . . ,  (A. 6) 
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10  6(uj2n-uj2n)  =  — 0.03f3  + 

+/[-7.20  sin(2A3  +  2D)  -  3.14  sin(2A3)  +  1.47  sin(A3  +  D  +  F)  - 
—  1.43  sin(2A3  +  2 D  +  /)  —  1.41  sin(A3  +  D  —  F)  +  0.29  sin(A3  +  D  +  F  +  /)  — 
-0.27  sin(2A3  +  4 D  -  l)  -  0.24  sin(2A3  +  4 D)  +  0.04  sin  3A3  +  0.18  cos  3A3  - 
-0.20  sin(2A3  +  2 D  +  21)  +  0.20  sin(2A3  +  2 D  -  l)  -  0.12  sin(A3  +  D  -  F  +  l)  - 
—0.11  sin(A3  +  D  —  F  —  l)  —  0.06  sin(2A3  +  4 D  +  /)  +  0.06  sin(A3  +  3 D  +  F  —  l)  + 
+0.05  sin(2A3  +  /)  +  0.05  sin(A3  +  3 D  +  F)  +  0.03  sin  2 (p  —  0.02  cos  2 (p  + 

+0.04  sin(A3  -\-  D  -\-  F  -\-  21)  —  0.04  sin(A3  +  D  +  +  —  /)]  + 

+t 2  [0.06  sin(2A3  +  2D)  -  3.51  cos(2A3  +  2D)  +  0.02  sin  2A3  -  1.53  cos  2A3  + 

+0.01  sin(2A3  +  2D  +  /)  -  0.70  cos(2A3  +  2D  +  /)  +  0.36  cos(A3  +  D  +  F)  - 
-0.34  cos(A3  +  D-F)-  0.13  cos(2A3  +  4  D  -  l)  -  0.12  cos(2A3  +  4  D)  - 
—0.10  cos(2A3  +  2D  +  21)  +  0.10  cos(2A3  +  2D  —  l)  —  0.08  sin3A3  +  0.01  cos3A3  + 
+0.07  cos(A3  +  +>  +  ++/)-  0.03  cos(2A3  +  AD  +  /)  -  0.03  cos(A3  +  +>-++/)- 
—0.03  cos(A3  +  D  —  F  —  l)  +  0.03  cos(2A3  +  /)  —  0.01  sin  A3  +  0.01  sin(A3  +  2D)  — 
-0.01  cos(2A3  +  2D  -  2  F)  +  0.01  cos(A3  +  3  D  +  F-l)]  + 

+73 [1.07  sin(2A3  +  2D)  +  0.08  cos(2A3  +  2D)  +  0.47  sin2A3  +  0.01  cos2A3  + 

+0.21  sin(2A3  +  2D  +  /)  -  0.09  sin(A3  +  D  +  F)  -  0.02  cos(A3  +  £>  +  +)  + 

+0.08  sin(A3  +  D  -  F)  +  0.04  sin(2A3  +  4D  -  l)  +  0.04  sin(2A3  +  4 D)  - 
—0.03  sin(2A3  +  2D  —  l)  +  0.03  sin(2A3  +  2D  +  21)  —  0.02  cos3A3  — 

-0.02  sin(A3  +  D  +  +  +/)]  + 

+74  [ — 0.04  sin(2A3  +  2D)  +  0.24  cos(2A3  +  2D)  +  0.11  cos  2A3  + 

+0.05  cos(2A3  +  2D  +  /)  —  0.01  cos(A3  +  D  +  F)  +  0.01  cos(A3  +  D  —  +)]  — 

—t5  [0.04  sin(2A3  +  2D)  +  0.01  cos(2A3  +  2D)  +  0.02  sin2A3]  +  ....  (A. 7) 


The  coefficients  of  the  series  (A.4)-(A.8)  are  expressed  in  radians. 


